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1 Flavor- Topology Duality 



Through the work of 't Hooft we know that the low-energy limit of QCD is dominated by 
light flavors and topology [|J . We expect that the same will be the case for the low- lying 
eigenvalues of the Dirac operator. Indeed, for massless flavors, the fermion determinant 
results in the repulsion of eigenvalues away from A = 0. It is perhaps less known that the 
presence of exactly zero eigenvalues has the same effect. The reason is the repulsion of 
eigenvalues which occurs in all interacting systems and has probably best been understood 
in the context of Random Matrix Theory where the eigenvalues obey the Wigner repulsion 
law §. 

In QCD, the fluctuations of the low-lying eigenvalues of the Dirac operator are de- 
scribed by chiral Random Matrix Theory (chRMT) 0, ^, [5]. This is a Random Matrix 
Theory with the global symmetries of the QCD partition function. It is characterized by 
the Dyson index || (3 which is equal to the number of independent variables per matrix 
element. For QCD with fundamental fermions we have (3=1 for N c = 2 and (3 = 2 
for N c > 2. For QCD with adjoint fermions and N c > 2 the Dirac matrix can be repre- 
sented in terms of self-dual quaternions with (3 = 4. The main ingredient of the chRMT 
partition function is the integration measure which includes the Vandermonde determi- 
nant. In terms of Dirac eigenvalues i\k it is given by Yik<i ~ ^f\ l3 ^k U+ ^^ 1 - Therefore, 
the presence of Nf massless flavors, with fermion determinant given by Yl k X k / , has the 
same effect on eigenvalue correlations as v = 2Nf/f3 zero eigenvalues. More precisely, the 
joint eigenvalue distribution only depends on the combination 2Nf + (3v ||, [7|]. Based 
on the conjecture p, Q that the zero-momentum part of the QCD partition function, 
ZN ft u{jnii • • • , mjv>), is a chiral Random Matrix Theory, it was suggested || that its mass 
dependence obeys the duality relation 

ZN s A m ii ' ' ' , m N f ) ~ Y[m' f 'Z Nf+ufi (m 1 , ■ ■ • ,m Nf ,0, • ■ ■ ,0). (1) 

/ 

This relation, which is now known as flavor topology-duality, is a trivial consequence of 
the flavor dependence of the chRMT joint eigenvalue distribution. 

The mass dependence of the chRMT partition function can be reduced to a unitary 
matrix integral which is known from the zero momentum limit of Chiral Perturbation 
Theory || W\. Starting from this representation of the low energy limit of the QCD 
partition function, also known as the finite volume partition function, flavor-topology 
duality was first proved for Nf = 2 in |J. However, its generalization to arbitrary Nf 
and v has only been achieved for (3 = 2 [II], [EJ. The relation ([]]) has been particularly 



useful for establishing relations between correlation functions of Dirac eigenvalues and 



finite volume partition functions [O, 13, 14, 15, 16 



For (3 = 2 the unitary matrix integral, which represents the low-energy limit of the 
QCD partition function, is also known as the one-link integral of two-dimensional QCD 
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(see for example fl7[l ) or the Brezin- Gross- Witten model [18|. For zero topological charge 
it can also be represented as a generalized Kontsevich model [19] with potential V(X) = 
1/X. This model has been discussed extensively in the context of topological gravity (for 
reviews see 



23] , |24| , |25|.) In this case the one-link integral has been analyzed |19| |2(| by 



means of Virasoro constraints which are based on its invariance properties. The main 
issue we wish to address in this article is whether flavor-topology duality of the one-link 
integral can be understood without relying on its chRMT representation. We will do 
this by deriving Virasoro constraints for arbitrary topological charge. For (3 = 2 we find 
a nonperturbative solution of these constraints in the form of a generalized Kontsevich 
model. 

In section 2 we will derive the Virasoro constraints for arbitrary v and three different 
values of (3. We observe that they satisfy the flavor-topology duality relations. A recursive 
solution of these relations is presented in section 3. It also provides us with an efficient 
derivation of sum-rules for the inverse Dirac eigenvalues to high order. In the second part 
of this section we discuss the uniqueness of the nonperturbative solution of the Virasoro 
constraints for (3 = 2. Concluding remarks are made in section 4. 



2 Virasoro Constraints 

In the simplest case, (3 = 2 and vanishing topological charge v — 0, the small-mass 
Virasoro constraints were first found in j0| after an identification of the appropriate 
unitary integral and the corresponding Generalized Kontsevich Model (GKM). Here we 
work entirely in the context of unitary integrals and derive a simple form of the Virasoro 
constraints valid for arbitrary topological charge and f3. They are obtained by expanding 
the partition function in powers of the masses. 



2.1 f3 = 2 

The low energy limit of QCD in the phase of broken chiral symmetry is a gas of weakly 
interacting Goldstone bosons. Its partition function is determined uniquely by the invari- 
ance properties of the Goldstone fields. For quark masses m and space-time volume V in 
the range 

I < v 1 ^ < — , (2) 

A 777™ 



(where m w = ^2mY,/F 2 is the mass of the Goldstone modes with F the pion decay 
constant and £ the chiral condensate) the partition function for the Goldstone modes 
factorizes |], 10] into a zero-momentum part, also known as the finite volume partition 



function, and a nonzero momentum part. For fundamental fermions and a gauge group 
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with N c > 3, the finite volume partition function in the sector of topological charge v is 
given by the following integral over the unitary group || 10 

Zt\MM ] ) = I dU (detuy e^ MU ' +M ' u \ (3) 

JU£U(N f ) 

where M. = MVE. Here and below we always take v > 0. The quantity M stands for 
the original unsealed quark mass matrix. The partition function is normalized such that 
Z?= 2 {M, Ptf) -> det v {M) for M -> 0. 

Especially for v = 0, the integral (|3]) has been studied extensively in the literature. In 
the context of lattice QCD it is known as the one-link integral [18, 21, HJ (see also [O 
for a review). Of particular interest is the fact that for Nf — > oo the partition function 



2£ =2 undergoes a phase transition [|Tj| from a small mass phase ( M. — > 0) to a large 
mass phase ( — ► oo ). The large mass expansion is asymptotic and its coefficients can 
be calculated, even for finite Nf, by expanding about a flavor-symmetric saddle-point. 
For Re(.M) — ■> oo such expansion can be carried out in a efficient way by means of the 
the (large-mass) Virasoro constraints found in P^| . For Im(A4) — > oo, at fixed value of 
Re(A / i), other flavor-nonsymmetric saddle points contribute to the partition function as 
well, and it is not known whether the expansion can be carried out by means of Virasoro 
constraints (for an explicit calculation including nonsymmetric saddle-points see |[26||). 
On the other hand, the small-mass expansion has a non-zero radius of convergence, and it 
can be determined by different methods for any finite Nf. For example, one can expand 
the exponential in (|3|) and calculate the corresponding unitary integrals systematically 



using a character expansion [27, 28]. Another efficient technique is again the use of the 



(small-mass) Virasoro constraints found in [19, |2(| |I2| . By generalizing such constraints 



to arbitrary topological charge v we will show that they naturally lead to flavor-topology 
duality. 

First of all, using the unitary invariance of the measure in ([5]) we can deduce the 
covariance properties of the partition function under a redefinition of the mass matrix 
M -> V- l M, 

{det M)- V Z^= 2 {M,M ] ) = {de%V- 1 M)- v Z^ 2 {V- l M,M^V). (4) 
This equation implies that (det M)~ V Z^ =2 (M, M ] ) is a symmetric function of the eigen- 



values of the hermitean matrix L a b = (M.M.') a b- Following |TJ| we can introduce the 
infinite set of variables 

fc-Lftp^I)*, (5) 

which are explicitly symmetric with respect to permutations of the eigenvalues and write 

(det M)-»Zt 2 (M,M^) = G u (t k ), (6) 
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where G v {tk) has the Taylor expansion 



G v (tk 



1 + aiti + a 2 t 2 + a n t 1 + 



(7) 



A simple consistency check on (||) is that unitary integrals are only nonvanishing if the 
powers of U and U* are equal. Following |19| we can find the coefficients of the Taylor 
expansion (|7|) from the differential equation 



d 2 Zt 



Sr 



cb g/3=2 



dMbadMlc 4 

which follows directly from @. AS a consequence, G v (tk) satisfies the equation 



(8) 



dMbadMac 

Using the chain rule 





82 KA-l 9 

r + vM ob 



dM 



G u {tk 



-G v {tk) 



(9) 



dMl 







E 

k>l L 

E 

k>l 



{MM'f^M' 
M^MM^f- 1 



4fc 



(10) 



- ab 



with dk = d/dtk we immediately obtain from (7) 



s=l 



cb 



where 



s-1 

^f =2 = E ^s-fc + E kt k9s+k + (Nf + U)d s , 8>1 

k=l fe>l 



(12) 



and our convention throughout this paper is that terms like the first one in the expression 
for C s vanish for s = 1. The operators C s obey a sub-algebra of the Virasoro algebra 



[C r ,C a ] = (r - s)C r 



(13) 



without central charge and r,s > 1. 

At fixed value Nf the are not independent, and the coefficients of the expansion (0) 
are not determined uniquely^. In (|TT|) the matrix elements of AiAi^ and the tk are not 
independent so that we cannot conclude from ([TT|) that the coefficients of (AiAi^) 3 vanish. 
Indeed, the corresponding equations for Nf = 1 and v — are inconsistent. However, 

1 Of course, they are determined uniquely up to the order that the tk are independent. 
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in order to fix the coefficients in (0) uniquely, we may supplement (|TT]) with additional 
equations. We do that by requiring that all coefficients of (M.M.^) S vanish, 

Ct 2 -6s,x]G v {t k ) = 0, s>l. (14) 

This procedure is justified provided that the all equations are consistent. Indeed, for 
Nf — > oo the matrix elements and the tf. are independent so that ([14]) must be valid. 
By inserting the Taylor expansion (0) into (|14]) we obtain an inhomogeneous set linear 
equations for the coefficients which depend on the parameter Nf + v. Inconsistencies can 
only arise for isolated values of Nf + v for which the homogeneous part of the equations 
becomes linearly dependent. This is indeed what happens in the example Nf = 1 and 
i/ = 0. 

Because of the commutation relations the constraints for s > 3 follow from the 
constraints for s — 1 and s — 2. The coefficients of the Taylor expansion (0) can be found 
recursively by solving the constraints for s = 1 and s = 2. This will be carried out in 
section 3 after deriving the Virasoro constraints for (3 = 4 and (3=1. 

2.2 (3 = 4 

Similarly, we now derive the Virasoro constraints for adjoint fermions ((3 = 4). For 
fermions in the adjoint representation, the zero-momentum Goldstone modes belong to 
the coset space SU(Nf)/ SO(Nf). These Goldstone fields are conveniently parameterized 
by UU T with U G SU(Nf). In the sector of topological charge u, we thus find the finite 



volume partition function |T0|, |29 



Z§=\M s ,MU = dU (detU) 2V e^ M ^ uu ^ +M s uut ) , (15) 

JueU(N f ) 

where V = N c v and M. s = MHV . In this case the mass matrix Ais is an arbitrary 
symmetric complex matrix. From Ai,s — ► V l AisV we obtain the transformation law 

(det M s r P Z§ = \M s , Ml) = (det V'MsV)^ ^{V'MsV, V^M^V*)- 1 ) . (16) 



This implies that {det Ais) v ' Zjt 4 (A4g, Ai s ) is a symmetric function of the eigenvalues 
sMl 



of AtsAi\. We thus have that 



(det.M 5 )-^= 4 (.M 5 ,.A4) = Gt\4) = 1 + oftf + a s 2 t s 2 + af x (tf) 2 + • • • , (17) 
where in analogy with (Rf) we have defined 



^M^l , as) 
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and G^j 4 (tf) will be determined by the differential equation 

Ft 2 7 /3=4 

i = 5 bc Z§= A , (19) 
dM Sba dMl ac 

After substituting the expansion (|T7D we obtain a differential equation for G^ -4 (tf) which 
can be written in the form of the Virasoro constraints, 

y/3=4f+S 



£f 4 -25 Sil Gt\4) = 0, 8 >1, (20) 



where 



t*T 4 = 2 E d%df-k + E + + ^ + *)&. , (21) 

fc=i fc>i 

and <9f = d/dt%. One can easily check that £f =4 and £f =2 satisfy the same algebra. 
2.3 0=1 

For QCD with fundamental fermions and N c = 2 the Lagrangian can be written in 
terms of fermion multiplets of length 2Nf containing Nf quarks and Nf anti-quarks. 
The chiral symmetry is thus enhanced to SU(2Nf). Since the fermion condensate is anti- 
symmetric in this enlarged flavor space, the coset space of the Goldstone modes is given by 
SU{2N f )/Sp(2N f ). This coset manifold is parameterized by UIU* where U G SU(2N f ) 
and I is the 2Nf x 2Nf antisymmetric unit matrix 



1 
-1 



(22) 



The partition function in the sector of fixed topological charge v is then given by inte- 
grating over the group manifold U(2Nf) instead of SU(2Nf), 



Z^ =1 (M A ,M A ) = dU (detU) v ei*l M * uw ' +M * uw *> , >, (23) 

JUeU{2N f ) 

where the mass matrix A4a = MVE is an arbitrary anti-symmetric complex matrix. In 
addition to the usual mass-term given by 



M 
-M 



(24) 



it contains di-quark source terms in its diagonal blocks. Below, in the calculation of the 
mass dependence of the partition function, the di-quark source terms will be put equal to 
zero. The covariance properties of Z^ =1 are given by 

(det MaY v/2 Z t = \M a, M A ) = (det V 1 M A V)- U /2 Z? =1 (V 1 M A V,V- l M a V 1 - 1 ). (25) 
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We thus have that (det SAa) V ^ 2 Z^ 1 (A4a, -Ma) is a symmetric function of the eigenval- 
ues of M. A S/V A . This results in the expansion 

(det^r/ 2 ^ 1 ^,^) = Gl=\t£) = 1 + aft* + 44 + atMf + ••-, (26) 

where 

The factor | in the above definition of takes into account that Ai is a 2Nf x 2Nf 



matrix for (3 = 1. Substituting (j26|) in the differential equation 



^Zt 1 , (2f 



we deduce the differential equations for G^ =1 {t k ) in the form of Virasoro constraints: 

(£?= 1 - = 0, (29) 



with 



£f 1 = \ E + E kt A k df +k + {N f+ u -- s -)dj , (30) 

z fc=l fc>l z z 

and d k = d/dt k . The £f =1 satisfy the same algebra as the two other values of (3 discussed 
in previous sections. 

3 Solving the Constraints 

In the first part of this section we discuss the recursive solution of the Virasoro constraints. 
Flavor-topology duality and nonperturbative solutions are discussed in the second part of 
this section. 

3.1 Recursive Solution of Virasoro Constraints 

Remarkably, in all three cases, (3 = 1, (3 = 2 and (3 = 4, the constraints ([14]) , (|20"D and 
( p5D can be written in the unified form: 

(£f-I<5 S)1 )G( a , T ) = 0, (31) 
7 

where (see (p)), (P) and @) 

=-£ dfe<9 s -fc + E + - [a + 7(2 - s) + s - 1] 9 a , (32) 

7 fc=i fe>i 7 



8 



and, motivated by ChRMT results of 



31 



, we have introduced the notation 



a = (JV>- 2)2/0 + + 1 , 

7 = VP- (33) 

The tk are defined as in (|5|), (|1|) and (^7]) for (3 = 2, 4, 1, respectively. From (|3T|) 
we conclude that in all three cases the Virasoro constraints only depend on Nf and v 
through the combination Nf + (3u/2 (v — > V for (3 = 4). Proceeding further, we can 
recursively solve ( |32| ) by substituting the Taylor expansion of the form (|7|) for G(a,7) in 
(|31"|). By treating the as independent variables, we find the following relations between 
the expansion coefficients: 

oo 

7 ^2 H n k+i + l)am-n k -in k+1 +in k+2 - + (a + j)(ni + l)a ni+ i n2 ... = a nin2 ... , (34) 
fc=i 

for s = 1 and 

oo 

(ni + 2)(ni + l)a ni+2 „ 2 ... + 7 H n k+2 + l)a„ 1 ...„ fc _i rifc+1 n fc+2 +i- 

fc=i 

+ (a + l)(n 2 + l)a nin2+ i n3 ... = 0, (35) 

for s = 2. (The subscript niri2 ■ ■ • of the coefficients a n in 2 ... is a shorthand for the parti- 
tion l ni 2™ 2 • • •.) All higher order Virasoro constraints are satisfied trivially through the 
Virasoro algebra. If we denote the level of the coefficients by n = J2k and and the 
total number of partitions of n by p(n), the total number of unknown coefficients at level 
n + 2 is equal to p(n + 2) whereas the total number of inhomogeneous equation (for s = 1) 
is equal p(n + 1) and the total number of homogenous equations (for s = 2) is equal to 
p(n). The total number of partitions satisfies the recursion relation f3"0|| 

p{n + 2) = p{n + 1) + p{n) — p{n — 3) — p{n — 5) + • • • 

- (-l)V„ + 2-^*)-(-l)Vn + 2-^t*) + .... (36) 

Since the number of partitions is a monotonic function of n, we have that 

pin + 2) <p(n + l)+p(n), (37) 

and the number of equations is always larger than or equal to the number of coefficients. 
In Table 1 we give the total number of coefficients and the total number of equations up 
to level n = 10. 

From normalization condition % = 1 we conclude that a\ = l/(« + 7). At level n = 2 
we obtain one equation each from ( [54] ) and (|3"5]), respectively, 

(a + l)a 2 + 2an = 0, 
7a 2 + 2(a + 7)an = a x — l/(a + 7) , (38) 
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Table 1: The total number of unknown coefficients (P(n)) and the total number of equa- 
tions (P(n — 1) + P(n — 2)) at level n. 



n 


P(n) 


P(n- 1) + P(n-2) 


n 


P(n) 


P(n - 1) + P(n - 2) 


1 


1 


1 


6 


11 


12 


2 


2 


2 


7 


15 


18 


3 


3 


3 


8 


22 


26 


4 


5 


5 


9 


30 


37 


5 


7 


8 


10 


42 


52 



which are solved by 

an 



(1 + a) 



2a(a + 7)(a + 7 + l) ' 



-1 . . 

a 2 = — — r . (39) 

a(a + 7)(a + 7 + 1) 



Up to order (MM^) 4 we find 
Z^ Nf (M,M^) = (det MT x 

l | Tr(MMt) (a + lXTr.M.Mt) 2 Tr(A^A^t) 



4(a + 7) 32«(a + 7)(a + 7 + 1) 32a(a + 7) (a + 7 + 1) 

2Tr(.M.Mt)3 (a - 7 + 2)Tr(MM t )Tr(AlM t ) i 



4 3 3a(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) 4 3 2a(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) 
[(a-7)(a + 3)+2] [Tr(.M.Mt)l 3 



4 3 6a(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) 

(7 - 5a - 6)Tr(MMt) 4 

4 5 a(a + l)(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) (a + 7 + 3)(a - 27) 

[2a 2 - 4a( 7 - 2) - 7 7 + 6] Tr(A4A / f t )Tr(.M.Mt) 3 
'4 4 3a(a + l)(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) (a + 7 + 3) (a - 2 7 ) 

[a 2 + a(5 - 37) + 2 7 2 - 7 + 6] fTr^A^t) 5 



4 5 2a(a + l)(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) (a + 7 + 3) (a - 2j) 

[3a 2 ( 7 - 2) - a 3 - 6 - 4 7 2 + 13 7 + a(-ll + 14 7 - 2 7 2 )] (Tr.M.M t ) 2 Tr(.M.M t ) 2 

4 5 a(a + l)(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2)(a + 7 + 3)(a - 27) 
[6 + a 4 + a 3 (7 - 37) - 257 + I87 2 + a 2 (17 - 2I7 + 2 7 2 ) + a (17 - 437 + 147 2 )] (TiMM^f 
H 4 5 6a(a + l)(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) (a + 7 + 3) (a - 2 7 ) 

+ C((7W7W t ) 6 )] . (40) 

For (3 = 1 the mass matrix Ai has been expressed in terms of the standard Nf x A7 mass 



matrix that occurs in the QCD partition function (see eq. (|24j)). For (3 = 2 and (3 = 4 
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the matrix M. is a complex Nf x Nf matrix as well, but for (3 = 4 it is symmetric. We 
have checked that all coefficients in (^) agree with the corresponding expansion for (3 = 2 
obtained in p0| . When the masses are degenerate (A4Ai^) a b = H 2 $ab exact expressions 
for Zl N are known for all three values of the Dyson index. For instance, in the special 
case Nf = 2 and (3=1 the finite volume partition function is given by |29] 

z ^ = 2 v (ll) ( " + " Wy+W m 

"' 2 fcSV 4 / k\(v + k)\(v + k + 2)'.(2v + k)\' 1 ' 

which is a convergent series for Z^ 1 . The reader can check that ([KJ) only differs from 
(p|) by the overall normalization factor 2 x ~ 2v j{v\{y + 2)!). We have also made similar 
checks for /3 — 4. The poles that appear in the expansion (f40|) at finite iVy, the so called 
de Wit-'t Hooft poles P2| , cancel after rewriting the expansion in powers of products of 
the eigenvalues of M^Ai. They can be regulated by choosing a noninteger value for the 
topological charge 0, |33|] . 

The expansion in terms of traces of the quark mass matrix can be directly compared 
with (n^2idet(D + m a )) Yang-Miiis from which we derive spectral sum- rules for inverse 



powers of the eigenvalues of the Dirac operator [ID|. With related to the eigenvalues 
±iAfc of the Dirac operator through = VEXk, the first six sum rules for m a = are 
given by 

1 




4(ct + 7) ' 

1 



16a(a + 7)(a + 7 + 1) 
a + 1 

16a(a + 7)(a + 7 + 1) 
1 



32a(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) 

a - 7 + 2 
64a(a 2 - 7 2 )(a + 7 + l)(a + 7 + 2) 

(a- 7 )(a + 3) + 2 



• (42) 

64a(a 2 -7 2 )(a + 7 + l)(a + 7 + 2) v ; 

The sum is over the positive eigenvalues only. It is clear that all sum-rules will depend on 
the topological charge through the combination Nf + (3u/2. The first formula in ( f£2"| ) had 
already been obtained in |29j from the low-energy QCD finite volume partition functions. 
For (3 = 2 our sum-rules agree with the ones derived in |2l| obtained by means of Virasoro 
constraints and flavor-topology duality. The remaining sum-rules have never been derived 
before from the finite volume partition functions for (3 = 1,4, but all except the fourth 
one have been obtained from ChRMT by means of Selberg integrals for arbitrary (3 |31 
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3.2 Flavor-Topology Duality 

So far for the sum-rules, we now return to the flavor-topology duality relations which is 
our main result. To this end we construct a [Nf + (3v/2) x (Nf + (3u/2) matrix Ai from 
the Nf x Nf original mass-matrix Ai by adding zeros. We thus have that Tr(.M.M^) fc = 
Tr(.M A4. ) k . The flavor-topology duality relation then follows from the observation that 
all coefficients in the small-mass expansion ( pECp depend on the number of flavors explicitly 
through the combination Nf + (3v/2 and can be written as 

Z? Nf (M,M"<) = (det M)» Z% >Nf+M2 (M,M), (43) 

where v — > V for (3 = 4. For the cases where (3v/2 is not an integer, the flavor-topology 
duality should be understood as an analytical continuation^ in v. Strictly speaking, we 
have only shown the duality relation ( fi3|) for the recursion relations for the coefficients 
of the small-mass expansion of the partition functions. To complete our proof we have 
to show that the Virasoro constraints uniquely determine all coefficients and that the 
expansion is a convergent series. The convergence of the small mass expansion follows 
immediately from the compactness of the unitary group. Indeed, for (3 = 2, the small-mass 
expansion can be resummed to the r-function of a KP-hierarchy |19|, . This result can 
also be obtained from an exact reconstruction of Z[j =2 [E?| from the character expansion 
2"7|| which is an alternative to the small- mass expansion derived here. The uniqueness of 



the solution of the Virasoro constraints is a much more difficult question ||34|| . We have 
checked explicitly to the order (.M-M^) 8 that this is indeed the case for all values of the 
parameters. 



For (3 = 2, however, the uniqueness follows from the explicit solution [[35[ of the 
equation 



d 2 N 



-Zt 2 (M,M^) = -±Z^= 2 (M,M X ) , (44) 



dM ba dM\ b J 4 

obtained from the trace of (H) or the Schwinger-Dyson equation ET 



r) 2 N 
M « M » *ka W Z^{MM) = MlM ab ^Zt 2 (M,Ml) , (45) 

OMbaOMhc 4 

obtained by contracting (|8|) with Ai and Ai^ . In both cases, the invariance of the partition 
function can be used to reduce the equations to a separable differential equation in the 
the eigenvalues of Ai^Ai. Because we are looking for a solution that is symmetric in the 
eigenvalues satisfying the boundary condition G u (0) = 1, it is determined uniquely by the 
solution of the single particle equation which happens to be the Bessel equation. 



2 The analytical continuation in v is subtle. The correct continuation is obtained by |26, [12| express- 
ing the finite volume partition function in terms of modified Bessel functions and make an analytical 
continuation in the index of the Bessel functions. 
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Let us discuss this in more detail for (|44[) . The first step is to write the partition 
function as 



detM_ 
detA^t 



u/2 



Z V {A) 



(46) 



where Z V (A) is a function of the eigenvalues Xf. of Jvft M. only (A = diag(xi, • • -x^ f )). It 
satisfies the differential equation 



y— 



dMtadMl 

This equation can be further simplified by introducing the reduced partition function 

z„(A) = A(4)nV^(A), (4 



(47) 



n 



k<l 



X i 



■1 j 



We thus 



where the Vandermonde determinant is defined by A(x 2 
have that z u (A) is a completely antisymmetric function of the eigenvalues with boundary 
condition such that Z U (A) ~ Yik x k f° r x k 0. The reduced partition function satisfies 
the separable differential equation 

4z/ 2 - 1" 



E 



si 



Ax 2 



z v {A) = N f z v {A) 



The solution can thus be expressed as a Slater-determinant 

z v {A) = det k i[(f) k (xi)], 



(49) 



(50) 



where the (j) k are the regular solutions of the single particle equation up to terms that 
vanish in the determinant. This requires further discussion. The solution of the single 
particle equation 

4z/ 2 - 1" 



D B 4>k{x) 
is the modified Bessel function 



0. 



Ax 2 



(j) k (x) = LO k (j)k(x) 



(f) k (x) = yfxl u (u k x). 



(51) 



(52) 



In order to cancel the Vandermonde determinant for x k — > we necessarily have that all 
uj\ are equal and thus u k = 1. Because of the antisymmetry of the columns in the Slater 
determinant (15(3) we can allow for single particle solutions that satisfy the equations 



D B <pi(x) = 4>i(x), 

D B (j) 2 {x) = <f) 2 (x) + /i 2 i0i(x), 

D B fa{x) = + HsKpiix) + H32<p2(x) ,etc. . 



(53) 
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By redefining 3 , 4 , • • • we obtain the equations (for k > 2) 

D B (f) k (x) = <p k (x) + n kk -i4> k -i(x) , (54) 

where the fj,kk-i<Pk-i(x) cancel in the expression for the Slater determinant. The solution 
of these equations follows immediately from the recursion relation 

D B [x k+ ^I u+k {x)} = x k+ ^I u+k (x) + 2kx k ~ 1 *I u+k _ 1 (x). (55) 

We thus find that 

fc = x k+ h v+k {x). (56) 

and /ifcfc-i = 2k. Since the Bessel equation has only one regular solution, we conclude that 
this solution is unique. The Virasoro constraints together with the boundary conditions 
determine the partition function uniquely and is given by a Slater-determinant of modified 
Bessel functions ([5(]) which agrees with the known result for the one-link intergral [ pl| , 

Actually the solution of the Virasoro constraints, although not its uniqueness can be 
obtained in a simpler way by identifying them with constraints which are satisfied by the 
Generalized Kontsevich Model (GKM). By changing variables from the matrix elements 
of Ai and Ai' to the matrix elements of the hermitean matrix L = (AiAi^/A) we can 
rewrite @ as 



d _ d d 
-L d e— h v- 



G V {L) = 6 bc G u (L) . (57) 



This equation admits a solution in the form of a GKM with potential V(X) = 1/X where 
X is a Nf x Nf normal matrix. More explicitly we haveQ, 



which is a simple modification of the v = case studied in ||19|| . In order to verify that 
([Blf) is a solution of (^) we have to insert a total derivative in the integral over X and 
multiply it by a second order differential operator as follows 

9 9 '"X^i^l^e^^^o. (59) 



dL k dL dc J aX ci \{AetX)"i 

Concluding, for (3 = 2 the Virasoro constraints naturally lead to a simple and rigorous 
proof of flavor-topology duality (|43D which is an alternative to the proof presented in 



3 The contour of integration is such that the integral converges for all values of L, and that an integral 
over a total derivative of X vanishes. 
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TT| , pL2| . For (3=1 and (3 = 4 a direct solution of the Virasoro constraints is not 
known, although we believe that further progress can be achieved at least in the case of 
equal quark masses where exact formulas for the partition function are known p9 |. For 
non-degenerate masses the obstacle is the lack of a generalization of the Itzykson-Zuber 
formula for integrals over unitary matrices. Therefore our proof of flavor-topology duality 
for (3 = 1 and (3 = 4 is only valid under the assumption of the uniqueness of the solution 
of the Virasoro constraints. 



4 Conclusions 

In a parameter range where the main contribution to the QCD partition function comes 
from the zero-momentum component of the Goldstone modes, the partition functions 
reduce to zero-dimensional integrals over the group manifolds SU(Nf) , SU(Nf)/0(Nf) 
and SU(2Nf) / Sp(2Nf) , for fermions in the fundamental representation and N c > 3 {(3 = 
2), N c = 2 (/3 = 1), and adjoint fermions with N c > 2 {(3 = 4), respectively. In the sector of 
fixed topological charge v we have shown that these partition functions satisfy differential 
equations in the form of Virasoro constraints which determine the small-mass expansion 
for arbitrary values of the topological charge. We have calculated the coefficients of the 
small mass expansion recursively up to the order m 8 and have obtained new spectral sum 
rules for inverse powers of the eigenvalues of the Dirac operator for (3 = 1, (3 = 2 and (3 = 4. 
For all three values of (3 the constraints depend on the topological charge through the 
combination Nf + (3u/2 which demonstrates flavor-topology duality for the corresponding 
small-mass (m — > 0) expansions. For the special case of (3 = 2 the uniqueness of the 
solution of the Virasoro constraints and the convergence of the small mass expansion 
implies flavor-topology duality for the full (finite m) partition function. In this case, 
we have found a non-perturbative solution of the Virasoro constraints in the form of 
a Generalized Kontsevich Model for arbitrary topological charge v which complements 
the v = results of We have not been able to find an analogous solution for 

(3=1 and (3 = 4. Based on numerical evidence we believe that for these two cases the 
solution of the Virasoro constraints is unique as well, but a rigorous proof could not be 
given. The existence of exact expressions for partition functions with degenerate masses 
{M.M^)ab = H 2 bab indicates that further progress is possible. 

It is worth commenting that, even for (3 = 2, although order by order the character ex- 
pansion is equivalent to the perturbative solution of the small-mass Virasoro constraints, 
it does not provide us with a natural proof of flavor-topology duality (see 123). It should 



also be pointed out that the GKM are known to satisfy both small-mass and large-mass 
Virasoro constraints. Thus, in view of our results, we also expect to find large-mass con- 
straints for v ^ at least for (3 = 2. Indeed we have explicitly obtained such constraints 



15 



entirely within the context of unitary integrals and they depend quadratically on the topo- 
logical charge v, contrary to the linear dependence of the small-mass constraints. The 
generalization of the large-mass constraints to (3 = 1 and (3 = 4 seems to be within reach. 
Due to the absence of an explicit Nf dependence of the expansion coefficients they are 



the most relevant ones from the point of view of both the replica limit |20], ^6|, [12] and the 
classification of universality classes . In particular, for applications in 2D-gravity they 
correspond to the continuum Virasoro constraints which appear after the double scaling 
limit. 
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